Calculations of collisions between cold alkaline earth atoms in a weak laser field 
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We calculate the light-induced collisional loss of laser-cooled and trapped magnesium atoms for 
detunings up to 50 atomic linewidths to the red of the ^So-^Pi cooling transition. We evaluate loss 
rate coefficients due to both radiative and nonradiative state-changing mechanisms for temperatures 
at and below the Doppler cooling temperature. We solve the Schrodinger equation with a complex 
potential to represent spontaneous decay, but also give analytic models for various limits. Vibrational 
structure due to molecular photoassociation is present in the trap loss spectrum. Relatively broad 
structure due to absorption to the Mg2^Eu state occurs for detunings larger than about 10 atomic 
linewidths. Much sharper structure, especially evident at low temperature, occurs even at smaller 
detunings due to of Mg2^ng absorption, which is weakly allowed due to relativistic retardation 
corrections to the forbidden dipole transition strength. We also perform model studies for the other 
alkaline earth species Ca, Sr, and Ba and for Yb, and find similar qualitative behavior as for Mg. 
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I. INTRODUCTION 



A. Background 



Laser cooling and trapping of neutral atoms has re- 
cently opened many new research areas in atomic physics. 
One can cool a gas of neutral atoms in magneto-optical 
traps (MOT) down to temperatures of 1 niK and be- 
low, and obtain densities up to 10^^ atoms/cm^. Evap- 
orative cooling methods have allowed the cooling of al- 
kali species to much lower temperatures below 1 /xK so 
that Bose-Einstein condensation (BEG) occurs. Binary 
atomic collisions play an important role in the physics 
of cold trapped atomic gases, and have been widely in- 
vestigated One of the first cold collisional process 
to be studied is the heating and loss of trapped atoms 
which result from tuning laser light to near-resonance 
with the atomic cooling transition Q . Here we take near- 
resonance, or small-detuning, to mean detuning A up to 
50 natural linewidths to the red of atomic resonance. 

Studies of small-detuning trap loss, extensively re- 
viewed by Weiner et al. , have mainly concentrated on 
alkali atoms ||] , for which it has been very difficult to de- 
velop quantitative theoretical models to compare with ex- 
periment. This is because alkali atoms have extensive hy- 
perfine structure, and thus the number of collision chan- 
nels is simply too large for accurate theoretical modeling. 
It has even been difhcult to estimate the relative weight 
of the different possible loss processes. Although one 
can develop simplified models, these are difficult to test 



adequately with complex alkali systems. On the other 
hand, trap loss photoassociation spectra in alkali systems 
for large detuning can be modeled quite accurately ■ 
This is because the underlying molecular physics of alkali 
dimer molecules is well-known, and the spectra are deter- 
mined by isolated molecular vibrational-rotational levels, 
for which the photoassociation line shapes can be well- 
characterized even in the presence of hyperfine structure. 
Quantitative analysis of such spectra have permitted the 
determination of scattering lengths for ground state colli- 
sions . These scattering lengths are critical parameters 
for BEG studies. 

Alkaline earth cooling and trapping have recently been 
of considerable experimental interest. Trap loss collisions 
have been studied in a Sr MOT |7| , and intercombination 
line cooling of Sr has resulted in temperatures below 1 /j,K 
and raised the prospects of BEG of Sr p|-pT||. Ga is of 
interest for possible applications as an optical frequency 
standard ||l2[-[l6| , and photoassociation spectroscopy in a 
Ga MOT has been reported Intercombination line 
cooling has also been reported for Yb [^,^, which we 
have included in our discussion because of its similarity 
in structure to alkaline earth atoms. 

Alkaline earth species provide an excellent testing 
ground for cold collision theories, especially given the 
rapidly developing experimental interest in the subject. 
Since the main isotopes of alkaline earth atoms have no 
hyperfine structure, the number of collision channels be- 
comes low enough to allow theoretical calculations even 
in the small detuning trap loss regime. Gonsequently, 
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this paper presents theoretical predictions for small de- 
tuning trap loss spectra in cold and trapped Mg gas in the 
presence of near-resonant weak laser light tuned near the 
^So — *■ ^Pi atomic transition, and discusses the nature 
of similar processes for Ca, Sr, Ba, and Yb. This work 



extends our previous note on Mg trap loss |20| to other 
species and lower temperatures, and shows the relation 
between small detuning trap loss and photoassociation 
theory. It is necessary to include spontaneous emission 
in modeling trap loss collision dynamics because of the 
long time scale of cold collisions. In addition, relativis- 
tic retardation effects play a prominent role at small de- 
tunings by allowing transitions to the dipole forbidden 
"'^Ilg state, which exhibits resolved vibrational structure, 
especially at very low temperature. Although we treat 
the long-range molecular interactions accurately, the po- 
tential energy curves and coupling matrix elements for 
the dimer molecules in the short-range region of chemical 
bonding are not sufficiently well-known to determine all 
aspects of trap loss. Therefore, we examine the uncer- 
tainties associated with the unknown phases developed 
in the short-range region of chemical bonding, and show 
which features are robust with respect to such uncertain- 
ties and which must be measured or later determined 
from improved theory. 



B. Trap loss collisions 

Light-induced trap loss takes place as a molecular pro- 
cess. Two colliding cold atoms form a quasimolecule, 
and their motion can be described in terms of the elec- 
tronic (Born-Oppenheimer) potentials of the molecular 
dimer with light-induced transitions between the molec- 
ular states. We consider only red detuning, which excites 
attractive potentials at long range. Such potentials sup- 
port a number of bound vibrational states. 

Figure |l| schematically indicates the nature of the trap 
loss process in a weak radiation field. An excitation laser 
with energy hi^ is tuned a few atomic linewidths below 
the atomic transition energy hvQ. The ground \g) and 
excited |e) quasimolecular electronic states are thus cou- 
pled resonantly by the laser at some long-range Condon 
radius Rc, where the photon energy matches the differ- 
ence between the excited and ground potential energy 
curves. The excited state decays to a loss channel \p) 
due to interactions at short range. The fully quantum 
mechanical description in Ref. |2^] shows that the overall 
probability Ppg of a trap loss collision can be factored as 
follows: 



P — P J 



(1) 



Here Jgg represents an excitation-transfer probability 
which is proportional to the scattering flux that reaches 
the short range region near Rp due to long-range exci- 
tation near the Condon point followed by propagation 



on the excited state to the short-range region. Ppe rep- 
resents the probability of a transition from the excited 
state to the loss channel at short range during a single 
cycle of oscillation through the short-range region. 

Figure |l] indicates the qualitative behavior of the trans- 
fer function J eg versus detuning A. At very small detun- 
ing of a few atomic linewidths Jeg becomes very small if 
there is a high probability of spontaneous emission dur- 
ing transit from the outer to inner regions. At sufhciently 
large detuning Jgg exhibits resonance structure due to 
the bound vibrational levels in the excited state poten- 
tial, where the vibrational period is much shorter than 
the decay time . This is the domain of photoassocia- 
tion to individual vibrational levels. Section III C below 
will show simple analytic formulas for J^g that apply in 
these two limiting cases of small detuning with fast decay 
or isolated photoassociation lines. These formulas show 
how Jf,g can in turn be factored as 



J — J P 

^ eg — 'J ee^ eg: 



(2) 



where Peg represents the probability of excitation from 
the ground to excited state in the outer zone near the 
Condon point, and Jee represents an excited state trans- 
fer function between the long-range outer excitation zone 
and the short-range zone. The factorization in Eq. (|^) is 
also schematically indicated in Fig. ^ 
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Transfer Excitation 




p c 
FIG. 1. Schematic representation of trap loss collisions. 

Each of the factors in Eqs. and (^) can be af- 
fected by unknown phases associated with the short- 
range molecular physics of the dimer molecules. (1) Peg 
is sensitive to the asymptotic phase of the ground state 
wave function. (2) Vibrational resonance structure in 
J< 

features. (3) Pj 
in short-range curve crossing p roba bilit ies. These effects 
are discussed in detail in Sees. II D and HE. The overall 



eg is sensitive to the position and widths of vibrational 
(3) Ppe is sensitive to Stiickelberg oscillations 



effect of such sensitivities will depend on the temperature 
and the alkaline earth species. 

There are two possible inner zone loss processes char- 
acterized by Ppe'- the state change (SC) and radiative es- 
cape (RE) mechanisms, both represented schematically 
by the loss channel p in Fig. 0. In the SC process the 
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excited state couples to another molecular state near 
a short-range crossing point Rp, and population trans- 
fer between them is possible. The products of the col- 
lision emerge on a state that correlates asymptotically 
with other atomic states of lower energy, such as '^P-f ^S, 
thereby releasing a large amount of kinetic energy to the 
separating atoms. In RE the excited state can decay by 
spontaneous emission after the atoms have been acceler- 
ated towards each other on the excited state potential. 
The ground state atoms then separate with this gain in 
kinetic energy. If enough kinetic energy has been gained 
to exceed the trap depth, this is a loss process (early de- 
cay after insufficient acceleration only leads to radiative 
heating). Here Rp represents the distance at which the 
atoms have received sufficient acceleration to be lost from 
the trap. Any emission with R < Rp leads to trap loss. 

In this paper we use the fully quantum complex po- 
tential method of Ref. |^l|] and do not resort to semiclas- 
sical methods (with one exception), although semiclas- 
sical concepts are often useful for interpretation. Our 
quantum methods are fully capable of describing the 
s-wave limit and the quantum threshold properties of 
the extremely low collision energies near or below the 
critical temperature for BEG. We describe the sponta- 
neous emission processes with a complex potential, and 
solve the corresponding time-independent multichannel 
Schrodinger equation in the molecular electronic state 
basis. This takes any vibrational state structure into ac- 
count automatically without the need to calculate wave 
functions or Franck-Condon factors, but limits our study 
to weak laser fields only, where only a single photon is ex- 
changed with the field. Typical cooling lasers are strong 
and detuned only a few linewidths. By alternating with 
a weaker probe laser one can access the particular range 
of detuning and intensity that we study. Future studies 
are needed to address the effects of strong laser fields and 
the consequent revision of our results due to saturation 
and power broadening p3| , p4[ . 

We also include rotational states into our model, which 
for the molecular ground state correspond to the partial 
waves of a standard scattering problem (angular momen- 
tum quantum number /). The symmetry of spinless al- 
kaline earth dimers permits only even partial waves. It 
should be pointed out that for near-resonant light the 
Condon point is at relatively large distances. This means 
that although the collision energy is low, one needs to go 
to relatively large / before the ground state centrifugal 
barrier stops the quasimolecule from reaching Re- We 
correct here some mistakes that we discovered in the sum 
over partial waves in Ref. |^ . 

The probe laser can be tuned over a wide range from 
a few to many atomic linewidths. For sufficiently large 
detuning the rotational structure becomes sharp, and it 
should be possible to resolve the vibrational and rota- 
tional states. However, even at 27 GHz detuning, rota- 
tional features in the Ca2 ^Tig — > ^Yiu photoassociation 
spectrum are only partially resolved [Q. Photoassocia- 
tion studies can yield precise information on the molecu- 



lar potentials. Also, the photoassociation line shapes are 
sensitive to the near-threshold ground state wave func- 
tion, especially if it has nodes in the region swept by 
the detuning-depcndcnt Condon point. Analysis of such 
data would hopefully give a value for the s-wave scat- 
tering length for Mg or other alkaline earth species, and 
consequently determine whether a stable Bose-Einstein 
condensate is possible. 

In this paper we present the calculated estimates for 
trap loss rate coefficients in Mg at temperatures around 
and below the Doppler temperature. The contributions 
from different mechanisms and states are identified and 
compared. We also calculate predictions for other al- 
kaline earth atoms and Yb by combining the appropri- 
ate atomic properties with model molecular potentials. 
Sec. H presents in detail the atomic data, molecular po- 
tentials, and laser couplings used for Mg and other alka- 
line earth atoms. Sec. Ill describes our theoretical ap- 
proach. The results for Mg are given in Sec. and for 
the other atoms (Sr, Ca, Ba, Yb) in Sec. Finally we 



present some conclusions in Sec. VI 



II. THE MOLECULAR PHYSICS OF ALKALINE 
EARTH DIMER MOLECULES 

A. Alkaline earth atomic structure 

Table || gives the basic atomic data for the alkaline 
earth atoms and Yb, which has an electronic structure 
similar to the Group II elements. The major isotopes 
have no hyperfine structure, except for Be. Alkaline- 
earth atoms have a ^So ground state, and excited ^Pi 
and ^Pi states that are optically connected to the ground 
state, the latter weakly. Figure g sketches the energy lev- 
els of the Group II atoms. The ^'^P first excited states 
are the most important for laser cooling. The ^'^D states 
shift downwards as the atomic number increases. For 
Mg the ^''^D states are above the ^■'^P states. For Ca 
and Sr they are between the ^P and the ^P states, and 
for Ba the ^''^D states are below the -'^'^P states. In laser 
cooling one uses typically the ^Sq-^Pi transition, which 
is the situation studied in this paper. For Mg this re- 
quires a UV laser source, and for the heavier elements 
requires repumping to recycle atoms that decay to lower 
levels. The weak -'^Sq-'^Pi intercombination transition has 
a very narrow linewidth, and is within the optical range. 
Thus alkaline earth atoms are good candidates for opti- 
cal atomic clocks when cooled to low temperatures. For 
clock applications we need to understand their laser cool- 
ing properties, including the magnitude and nature of 
laser-induced colhsional trap loss. 

Beryllium is not likely to be a serious candidate for 
laser cooling. It is toxic, has a very short wavelength cool- 
ing transition, and the intercombination "clock" transi- 
tion is so weak as to be effectively forbidden. Therefore, 
we will not consider Be in the rest of this paper. 
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TABLE I. Atomic data for major isotopes of Group II elements and Yb without hyperfine structure (natural abundance 
shown). Most data is derived from [^. The lifetime r for 'Pi is taken to be the inverse of the 'Pi-'So spontaneous emission 
rate Fat/^i, thus neglecting weak transitions to other states for Ba and Yb. The '^P lifetimes are from several sources [Q. The 
linewidth in frequency units is Fat/ft = {2-!tt)~^ . The wavelengths A and fine structure splittings for Sr and Ba are taken from 
[ p7[ . The Doppler temperature is defined as Td ~ rat/(2fcs). We take the recoil temperature to be Th = {h'^)/{mX^kB), where 
A = A/27r and m is the atomic mass. The dipole moment is do ~ \/ SFat A^/4 (in a.u.). The atomic units for dipole moment, 
length, and energy are eao = 8.4783 x 10"^° Cm, ao = 0.0529177 nm, and eV(47reoao) = 4.3597482 x 10"'** J respectively. 





Be 


Mg 


Ca 


Sr 


Ba 


Yb 


Major isotopes without 
hyperfine structure 
(abundance) 


^Be(100%) 


^''Mg(78.99%) 
2''Mg(11.01%) 


*°Ca(96.94%) 


****Sr(82.58%) 
»«Sr(9.86%) 


'^**Ba(71.70%) 


'^*Yb(31.8%) 
'^2Yb(21.9%) 


r 

'Pi (ns) 
^Vx (ms) 


1.80 


2.02 
5.1 


4.59 
0.48 


4.98 
0.021 


8.40 
0.0014 


5.68 
0.00088 


Fat//l 

'Pi (MHz) 
^Pi (kHz) 


88.5 


78.8 
0.031 


34.7 
0.33 


32.0 
7.5 


18.9 
120 


28.0 
181 


Doppler-cooling limit 
'Pi (mK) 
^Pi (nK) 


2.1 


1.9 
0.75 


0.83 
8.0 


0.77 
179 


0.45 
2.8 10^ 


0.67 
4.4 10^ 


Recoil limit 
'Pi (^K) 
'Pi (/^K) 


39 


9.8 
3.8 


2.7 
1.1 


1.0 
0.46 


0.45 
0.22 


0.69 
0.36 


do (a.u.) 

'So-'Pi 


1.89 


2.38 


2.85 


3.11 


3.16 


2.35 


A (nm) 
'So-'Pi 


234.861 


285.21261 


422.6728 


460.733 


553.548 


398.799 


A = A/(27r) (ao) 
'So-'Pi 


706.4 


857.8 


1271.2 


1385.7 


1664.9 


1199.4 


FS splitting 
^P2-'Po (a.u.) 


1.36 10"^ 


2.77 10"* 


7.20 10"'' 


2.65 10"=^ 


5.69 10"^ 


1.10 10"^ 



Because of the lack of hyperfine structure the ba- 
sic laser cooling mechanism for alkaline earth atoms is 
Doppler cooling, for which the temperature limit Td, de- 
fined in the caption of Table |, is set by the linewidth 
Fat of the cooling transition (widths in this paper are ex- 
pressed in energy units, so that the decay rate is Fat/Ti). 
The lifetime of the alkaline earth -'^Pi state is between 
1.8 and 8.4 ns, giving a Doppler-cooling hmit between 
2.1 and 0.45 mK for the elements in Table |[ On the 
other hand the ^Pi state has a long lifetime with Doppler- 
cooling limits in the nK range. This can be compared to 
the photon recoil limit Tr, defined in the Table caption; 
Table | shows that is between 0.2 and 4 /iK. Thus the 
recoil limit is above the intercombination line Doppler 
cooling limit for Mg and Ca, nearly coincident with it for 
Sr, and below it for Ba and Yb. 

Although Sisyphus cooling and magnetic trapping is 
not available for ^Sq atoms as it is for alkali atoms, in- 
tercombination line cooling is possible for some Group 11 
species. This has been used to cool Sr to ^ 400 nK with 
relatively high phase space density > 0.1 If com- 

bined with far off-resonant optical traps and evaporative 
cooling it may become possible to obtain Bose-Einstein 
condensation with optical methods alone. 
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FIG. 2. Energy levels of Group II atoms and Yb. 
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B. Alkaline earth dimer molecular structure 

Figure || shows the lowest electronic potentials for 
Mg2 . There are only two states with attractive long- 
range potentials correlating with ^Pi + ^Sq that can be 



interatomic distance R with a magnitude on the order of 
twice the atomic linewidth. Thus the vibrational levels of 
the state near the -'^S+^P dissociation limit overlap 
strongly, and one does not expect to resolve them. One 
interesting point is that the dipole- forbidden ^Ilg — ^Y]g 



resonantly coupled to the ground state by laser light, transition becomes allowed at large^i? due to retardation 
namely and ^Ilg. Both states offer the possibil- 

ity for the SC and RE mechanisms. When comparing 
Figs. H and || one can see that Mg is special. For other 
alkaline earth atoms the molecular state picture is fur- 
ther complicated by the atomic D states below the ^Pi 
state. This increases the number of molecular states and 
thus the number of energetically available exit channels. 
The small number of molecular states is one reason we 
have chosen Mg as the basis for our studies. Theoretical 
calculations require precise information about the molec- 
ular potentials and couplings over a wide range of inter- 
atomic distance. Ab initio calculations of ground and 
excited molecular potentials are also available for Sr ||2^] 
and Ba (30|. 

While working on the manuscript we received new ab 
initio data on Mg2 from E. Czuchaj, including both im- 
proved potential curves and spin-orbit couplings pl| . Al- 
though the new data differ to some extent from the re- 
sults of Stevens and Krauss because of improved cal- 
culations of electron correlation effects, we do not expect 
that use of the new data would lead to any strong mod- 
ification of our basic results, which should be viewed as 
qualitative model calculations for the reasons to be dis- 
cussed in the following sections. 



corrections. This means that the ^Ilg state can be ex- 
cited at Rc, but the spontaneous emission probability 
goes down quickly as R decreases. Consequently, ^Ilg 
vibrational states near the dissociation limit have nar- 
row emission linewidths. Thus the assumption that the 
vibrational states overlap strongly and can not be re- 
solved at detunings of a few linewidths is not necessarily 
valid. One must determine if resolvable features persist 
when one sums over all involved rotational states/partial 
waves, and takes the energy average over a thermal dis- 
tribution. We show that one can indeed see vibrational 
structure, especially if the temperature is well below the 
^Sq-^Pi Doppler limit. As mentioned above, this is by 
no means impossible, if one does the cooling using the 
^So-'^Pi transition. 
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C. Long-range properties of states correlating with 
^Pi+^So atoms 

There are four molecular states correlating with 
^Pi-(-^So atoms: two long-range attractive states, 
and "'^Ilg, and two long-range repulsive states, and 
^n^. Here long-range means that R is large compared to 
the short-range region of chemical bonding and van der 
Waals interactions shown in Fig. || so that the potential 
is determined by the first-order dipole-dipole-interaction 
with C3 = T2rf§ and =Fc^o '^'^ ^^^^ ^11 states, re- 
spectively. Here do is the z-component of the atomic 
transition dipole matrix element, which is related to the 
atomic linewidth Fat and the reduced wavelength of the 
^So-^Pi transition (A = A/27r) by = 3A^rat/4. The 
exact long-range (Ir) potentials including relativistic re- 
tardation corrections are [02l 



STat 



[cos(m) + usin(u)], (3) 
[cos(w) -I- usin(u) — cos(m)]. 



FIG. 3. The molecular states of Mg2 in atomic units cor- 
responding to the asymptotic atomic states ^Sq-I-^Pi and 
^So4-^Po,i,2 [Q; the zero of energy is at the ground state 
^Sq-I-^So asymptote. There are four states correlating with 
each asymptote, of which two are attractive and two are re- 
pulsive at large K, where the system is expected to be resonant 
with the laser field. 

The linewidth of the excited ^Y^u state depends on the 



where u = RjX is the scaled distance. The molec- 
ular linewidths with relativistic retardation corrections 
are || 

Viu- ^ Fat |l - ^ [wcos(u) - sin(M)]l , (4) 



F(u; iHg) = Fat 1 1 - ^ {u cos(w) - (1 - u^) sin(u)] 



In the region with R < X, the potentials vary as 

and F(^S+) and F(-^Hg) respectively vary as 2Fat and 

Fat"V5- 
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The very long-range excited state potentials result 
in large Condon points for excitation of the attractive 
states. If the laser detuning relative to the atomic tran- 
sition is expressed in units of Fat and the distance in units 
of A, the scaled Condon point (uc) for the ground-excited 
state transition becomes independent of atomic species. 
Table |l| shows uc for several detunings A, where we de- 
fine red detuning to be positive. For Mg at A=rat we 
have i?c(^S+) = 1132 ao and Rc{^Ilg)^728 &o. 

TABLE II. Condon points in scaled distance for selected 
detunings 



Detuning A 






Fat 


1.32 


0.849 


5Fat 


0.716 


0.511 


lOFat 


0.555 


0.411 


30Fat 


0.376 


0.288 



The attractive potentials support a series of vibrational 
levels leading up to the dissociation limit. Assuming a 
potential with a long-range form —C^/R^ gives the bind- 
ing energy of vibrational level v p33[: 



.1/3 _ 



TT 

2^ 



[v - vof 



(5) 



where = (V¥/2)F(5/6)/F(4/3) = 1.120251, vd is the 
vibrational quantum number at the dissociation limit {vd 
is generally nonintegral) and is the reduced mass. The 
vibrational spacing function, which we need later, is 



dv 



— hv,, — 



as 



1/2 



-5/6 



(6) 



where hv^ is the vibrational frequency. 



D. Ground state 

The ground-state van der Waals potential varies at 
long-range as R~^ and is essentially flat for the range 
of Condon points we consider. The energy-normalized 
ground state scattering wave function for coUisional mo- 
mentum Tikoo and partial wave I has the long-range form 



^{R,l,k^) = 



2^ 



sm 



^R- -1 + rji{koo) 



(7) 



We define the coUisional energy as e = fi k^/{2ij). The 
short-range potentials are not sufficiently well-known for 
alkaline earth dimers to determine accurately the scatter- 
ing phase shifts ??;(fcoo). Therefore, our calculations will 
have to be model calculations. However, we test the sen- 
sitivity of our trap loss spectra to the unknown phases. 



and show that this is not a serious limitation. There are 
two reasons for this. One is that the ground state poten- 
tial is flat in the long-range region, and the amplitude of 
^(i?, Z,fcoo) has its asymptotic value in Eq. (M) indepen- 
dent of R as long as i? > = \{2^lCQ/h^Y/^ 0,11; 
the use of xq (or a closely related length) as an appropri- 
ate length scale for van der Waals potentials is described 
in the Appendix of Ref. |^ . The condition Rc ^ xq is 
easily satisfied in our case. Second, at the Doppler limit 
Td for ^So ^Pi cooling, a number of partial waves I 
contribute to trap loss i n our detuning range (1-50 Fat). 
We demonstrate in Sec. EVA that a sum over I removes 



the dependence on the short-range potential. 

Our trap loss spectra for s-wave scattering at the low 
temperatures available via intercombination line cooling 
will be sensitive to the actual scattering length Aq of the 
ground potential. However, we demonstrate a sim- 
ple scaling relationship that will allow our low T s-wave 
results to be scaled to any value of the scattering length. 
We may expect that the ^11+ scattering length can be 
determined from one or two-color photoassociation spec- 
tra, as has been done for alkali species However, such 
analysis will require an accurate Cg coefficient and for 
optimum results needs a reasonably accurate short-range 
potential as well. 



E. Excited state short-range potentials 

Trap loss spectra depend strongly on the excited state 
short-range potential structure in three ways: 

(1) The curve crossings leading to SC trap loss occur 
at short range, and determine Ppe- In a Landau-Zener 
interpretation. 



Ppe(J')=4e-2'^^ 



1 



)sin2(/3j0, 



(8) 



where A — \Vpe{Rp)\'^ / {TiVpDp) and J' design ates excited 
state total angular momentum (See Section II F). Here 
Vpe{Rp), Vp, Dp are the respective coupling matrix ele- 
ment, speed, and slope difference at the Rp crossing, and 
(3 J' is a semiclassical phase angle [Ba : 



(3j> = I k,iR,J')dR 



kp{R,J')dR + -, (9) 



where i?oi and hki{R,J') are the respective inner clas- 
sical turning point at zero energy and local momentum 
for state i = e or p. The attractive singlet potentials 
may have one or more crossings with repulsive potentials 
from lower-lying states, e.g., those correlating to ^P+^S. 
Thus the number, positions of crossings and the coupling 
between states with crossing potentials are important for 
the magnitude of Ppe ■ Although we can make reasonable 
estimates for A, the phase angle f3j' is sensitive to details 
of the potentials and can only be calculated accurately if 
very accurate potentials are available |M . 
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(2) The vibrational structure in the trap loss spectra 
depends on both short and long-range potentials. The 
spacing between the vibrational levels given by Eq. (||) 
depends only on the long-range potentials, but the exact 
positions of the levels in Eq. (j^) are determined by the 
short-range potentials in the region of chemical bonding 
(through the vd parameter). Thus, the magnitude of the 
vibrational spacings in our model calculations will be cor- 
rect, whereas the actual positions can only be determined 
by measurement. 

(3) The short-range SC process also contributes to the 
width of vibrational features in the trap loss spectra . 
Depending on species, temperature, and detuning, the 
widths may be primarily determined by natural or ther- 
mal broadening or by the predissociation decay rate re- 
lated to Ppe- However, we are able to plac e app roximate 



bounds on the magnitude of Ppe ■ Sections pi Q , |V|, and 
discuss the contributions to feature widths and show 
that natural and thermal broadening tends to be domi- 
nant at small detuning, whereas SC broadening may be- 
come dominant at large detuning. 

The available data on short-range potentials varies 
among the Group II elements. Ab initio potentials are 
available for Mg |^,^ . The structure of the molecular 
potentials is fairly simple because only states correlat- 
ing to ^P-I-^S (4 states) and ^P-|-^S atoms (4 states) are 
present (see Fig. ^) and the potentials provide quali- 
tative data for possible SC mechanisms. On the other 
hand, Ca, Sr and Ba have a very complex short-range 
structure, because of large fine structure splittings in the 
triplet states, and the states correlating to ^D-f^S and 
•^D-I-^S coming into play. For example, Boutassetta et 
al. |S] and AUouchc et al. 30 1 calculated the short-range 



potentials for Sr and Ba respectively, but the large num- 
ber of states (e.g., 38 states for Ba) makes it excessively 
complicated to treat the short-range SC mechanism even 
qualitatively. 

The coupling between the singlet and triplet states at 
short-range is due to spin-orbit couplings. The exact 
magnitude of the couplings is unknown for all Group 
II elements, but can be estimated using Table Al of 
Ref. ]39| , which relates the coupling matrix elements to 
the '^P2— '^Po fine structure splitting. This approxima- 
tion ignores any i?-dependence of the spin-orbit matrix 
elements due to chemical bonding. The fine structure 
splitting of the ^P state increases with atomic number 
(see Table I). 

We have opted for Mg as a model system, because 
short-range potentials are available. Each singlet state 
with attractive potential couples to only one triplet state 
in the inner zone. Thus the SC trap loss problem for 
Mg decouples into two three-state calculations, one for 
the and one for the ^Ilg excited state. For Mg we 
provide qualitatively correct SC trap loss spectra and 
even give some quantitative estimates. 

In the case of the other Group II elements and Yb, 
we treat their complicated inner zone physics as one ef- 
fective crossing, that is, we use three-state calculations 



based on the Mg-model, and explore the effects of mass, 
radiative properties, and coupling strength (size of spin- 
orbit splitting) in these model calculations. Thus we do 
not use the ab initio potentials of Refs. (2^j3^ discussed 
in the paragraph above, because even if we were to in- 
clude all the curves, there would still be uncertainties 
associated with unknown phases and unknown coupling 
matrix elements. Rather, our goal is to indicate qualita- 
tive differences in magnitude and spectral shapes among 
the various species. We trust that these will be helpful 
in providing guidance for future experimental and the- 
oretical studies of these systems. As we discuss in Sec- 
tion III B| , our calculation of the excitation-transfer coef- 
ficient K will be to a large extent independent of details of 
the complicated short-range molecular physics and curve 
crossings. 



F. Molecular rotational structure and coupling to 
the laser field 

The full three-dimensional treatment of the collision 
of a ^S atom coupled to a ^P atom by a light field is 
worked out in Ref. |Q. We adapt this treatment to our 
simplified model with molecular Hund's case (a) transi- 
tions between two molecular states with the usual rota- 
tional branch structure. The angular momentum J" in 
the ground state can only be that of molecular axis rota- 
tion, J" = I" with projection m" on a space-fixed axis. 
The excited rotational levels in a Hund's case (a) molec- 
ular basis for a or ^Hg state do not have mechanical 
rotation V as a good quantum number but instead the 
total angular momentum J' with space projection M' . 
The three possible transition branches have J' = I" + B, 
where the branch labels P, Q, and R respectively desig- 
nate the cases B = -1, 0, and 1. The quasi-molecule 
ground state can couple to the state only by P 

and R branches, but to the ^Ug state by P, Q, and R 
branches. All potentials have a centrifugal term added, 
n{n + l)/(2/ii?^), depending on ground state n — I" or 
excited states n — J' — I" + B. 

The radiative coupling terms are 



VA^B,l",m",qiu) 



_ /2ttI 



\ c 



1/2 



{A, J'M'\eq-d\l"m"), 



(10) 



where A — or ^Hg labels the excited state, / is light 

intensity, d is the lab frame dipole operator, eg is the 
polarization vector of light with polarization q = 0,±1, 
and M', m" are lab frame angular momentum projection 
quantum numbers. In our weak-field case, where transi- 
tion probabilities are linear in | V^^s, ;",,„"(«), gP) we can 
define an m"-averaged radiative coupling matrix element 
(the sum over m" removes the dependence on q): 

VegAu,l\B,I) 
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\ 7n — — I" 




= 2.669 X 10-^aA,B,i"\/l(W/cin)dAiu){a..u.). 

(11) 

The molecular electronic transition dipole moment cIa is 

(12) 



dA{u){a.u.) = / ^^(^■^■) rA(M)(a.u. 



where Ta{u) is the molecular linewidth of the excited 
state as in Eq. (Q). The factors \/2l" + 1 aA,B,i" are given 
in Table 

TABLE III. The rotational line strength factors 



V2l" + 1 aA,B,i- 



State A 



Branch B 



I" = (s-wave) 



E 
S 

n 
n 
n 



p 

R 

p 
Q 

R 










va"-i)/3 
y m" + 1)/3 

v/('" + 2)/3 



G. Model potentials for Mg 



Our results are not sensitive to the detailed form of the 



ground state potential, for the reasons given in Sec. [ID 
Thus, we model the ground state potential by a Lennard- 
Jones 6-12 form 



12 / cr 

XuJ \Xu 



(13) 



For Mg we model the ground state potential from ||2|] 
with a well depth of e = 0.002825 a.u. and an inner turn- 
ing point of (7 = 6.23 ap. This potential has a scat- 
tering length of —95 ao for the ^''Mg reduced mass of 
23.985042/2 atomic mass units. The scattering length 
Aq is determined from the fc ^ behavior of the s-wave 
phase shift: 770 = — fcoo^o- 

Since the excited ah initio potentials of Ref. ||2^ do not 
permit a quantitative calculation of s pectroscopic accu- 
racy, for the reasons given in Sec. [IE, the specific forms 
of the short-range (sr) excited-state potentials are not 
important for our purposes of modeling the qualitative 
structure and magnitude of the coUisional loss. How- 
ever, it is important to retain the correct long-range form. 
Therefore for simplicity in the calculations and because 
we want to model the other alkaline earth systems to ex- 
plore the effect of different C3, Fat, A and mass, we have 



modeled the ah initio potentials with Lennard- Jones 3-6 
potentials keeping the long-range form fixed to its known 
C3 value given in Eq. (^) 



\Xu) \\u) 



(14) 



We have two fitting parameters e and ct, and three 
given values: well depth of the ah initio potential (the 
depth of the model potential is e), position of minimum 
Mmin and C3(-^E+) — — 2dQ or C^{^Iig) = —d^. Because 
we want to fix the long-range potential C3 we can not 
fit Mmin and the well depth at the same time and have 
chosen the latter: 



K(zi;'S+)=4e(E) 



a(S) 
Am 

cos(m) 

a(n) 



(15) 



I sin(u)], 



"0 
"a\3 



Xu 

cos(m) -f usin(u) — ti^ cos(u)] 



(16) 



where 

e(S) = 0.0347 a.u., cr(S) = 4.339 ao, 

e(n) = 0.0681 a.u., cr(n) = 2.751 ao. 
The well minima are Awmin = 5.5 ao and 3.5 ao for 
and ^Ilg, respectively, compared to the ah initio values 
= 6.1 ao and 5.4 ao. 

Both excited states have a SC mechanism in the short- 
range region with coupling to a triplet state. The triplet 
states are purely repulsive, and modeled with 



X 



(17) 



where A — or ^Ilg labels the molecular state, and 
= -0.0601 a.u. 

The SC from to '^Hu takes place around the inner 
turning point of the ^T,^ potential well. We have chosen 
Cg = 392 a.u. for the model of the '^11,1 state. 

The SC from ^Ilg to '^E^ takes place about 1.5 ao out- 
side and 0.019 a.u. above the minimum of the ^Ilg state 
potential well. With Cg = 81 a.u. we have a model of the 
crossing where the corresponding values are: 1.0 ao and 
0.019 a.u. The difference in slope of the crossing poten- 
tials is 0.030 a.u. for the ah initio potentials and 0.037 
a.u. for the model. 

The coupling between the crossing states are approx- 
imated using Table Al of Ref. |39|. The ^S+-3n„ and 
^Ug-^T,'^ matrix elements are C/v^ and (/2, respectively, 
where C = 1-84 x lO""* a.u. is 2/3 of the atomic ^Pj- 
•^Po splitting. For example, we estimate an upper bound 
(sin^^j' = 1) to the Landau-Zener version of Ppe in 
Eq. (g) for the ^Ilg-^Sg crossing to be 2.6x10^^ for the 
model potentials and 2.8x10"'^ estimated from the ah 
initio potentials. This upper bound is consistent with 
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the calculated Ppe as a function of J' from our complex 
potential calculation described below. 

Since we will use a semiclassical method to determine 
the Ppe fa ctor for the RE process via the state (see 



Sec. 



HID 



below) , we do not need an explicit probe chan- 



nel for RE. However, we introduce a probe channel to 
simulate RE in order to show that the same Jeg factor 
in Eq. (0) apphes for both SC and RE processes, irre- 
spective of the choice of the short-range Rp. Since we 
may take any form we like for a RE probe state, we use 
a probe state potential which crosses the excited state 
at a distance Up, where the kinetic energy gained by the 
collision pair is 1 K, corresponding to a trap depth of 0.5 
K. The RE probe potential has a repulsive inner wall 



RE, probe 



C 



12 



[Xu) 



12 



in,RE,ctt 1 



(18) 



where Vkin,RE,oo = 3.17x 10~^ a.u. and C12 = 5x 10^ a.u. 
No rotational term is included in this probe channel. The 
same probe state potential is used for all collision ener- 
gies, which are small (mK range and below) compared to 
the 1 K kinetic energy at Up. The coupling between the 
excited state and the probe state is chosen to be weak: 



10- 



a.u. 



H. Model potentials for Ca, Sr, Ba and Yb 

Since the different ground state values of Cq and inner 
potential shape make no differen ce for these model stud- 
ies, for the reasons given in Sec. II G, we take the same 
ground state Lennard- Jones 6-12 potential, Eq. (|l3|), as 
in the Mg case to model the other alkaline earth ground 
states, and only change the reduced mass. This proce- 
dure yields respective model scattering lengths of 67 ao, 
-65 ao, -41 ao, and 97 ao for ^OCa, ^^Sr, ^^SBa, and ^^^yb. 
Thus, |v4o| <C Rc in all model cases. 

The long-range of the excited state potentials for Ca, 
Sr, Ba and Yb is still exact, using the data from Table 
I with the form in Eq. (^) . Due to lack of accurate ex- 
cited state short-range molecular potentials and because 
of their more complicated structure, we model the trap 
loss for Ca, Sr, Ba and Yb by scaling the potentials from 
the Mg model Eqs. (|l|) and (^. The well depth e of 
the ^Sj and ^11^ potentials are scaled by the size of the 
singlet-triplet states splitting compared to that splitting 
in Mg, e.g.: 



eca 



P(^P,Ca)-£:(3p,Ca) 
P(ip,Mg) - £;(3p,Mg)' 



(19) 



The short-range structure is treated as one effective 
crossing. The probe states are qualitatively like those in 
the Mg model. The position (in energy) of the crossing 
between the ^11^ and probe state potentials is scaled as 
above. The and probe state potentials come very 
close at the inner wall of the potential around the 
classical turning point K(u) — e. 



The spin-orbit coupling constant C scales with the spin- 
orbit splitting of the ^P atomic states. We use the same 
definition of the matrix elements as for Mg. The Landau- 
Zener adiabaticity parameter 27rA in Eq. (||) for Ba and 
Yb is larger than unity, leading to a modified shape of the 
short-range adiabatic potentials and very small Ppe,sc- 
Thus for Ba the ^E+-probe state coupling and for Yb 
the ^Ilg and the ^I]+-probe state couplings have been 
reduced by about a factor 5 to obtain values of Ppe,sc 
close to unity, in order to test the limit of very strong 
broadening of the vibrational structure. We believe this 
limit is physically more realistic. The variety of crossings 
in these systems might lead to a strong SC process. 



III. COMPLEX POTENTIAL CLOSE COUPLED 
CALCULATIONS 



A. Description of method 



The weak field approximation assumed in this study 
allows us to apply a complex potential method pl| , ^ , 
since re-excitation of any decayed quasimolecular popu- 
lation can be ignored. Furthermore, the weak field only 
couples each ground-state partial wave to at most 3 rota- 
tional states of the or ^11^ state through the P, Q, or 
R branches. However, in the weak field the excited-state 
rotational states do not couple further to other ground- 
state partial waves. Therefore we can ignore any partial 
wave ladder climbing. Thus, for each trap loss mecha- 
nism we have three dressed states: a ground state g, an 
excited state e, and a probe state, p. We solve the three- 
channel, time-independent radial Schrodinger equation 
for ground state collision energy e, partial wave I — I" , 
for each transition branch B and for a given intensity / 



(£,i?) 



2^ 



2 - V(i?, /, B, A, /)] 0(£, R) = 0, 



V is the 3x3 potential matrix 
V{u,l,B,A,I) = 



(20) 



(21) 



VegiuJ,B,I) 
Vpeiu) 



VgiuJ) 





Vp{u,l,B) 



The elements of V are described in Sec. ^ above. A 
complex term —iT(u)/2 is added to the excited-state po- 
tential to simulate the effect of excited-state decay. The 
full retarded form of the molecular linewidth, Eq. (jj), is 
used. 

Application of standard asymptotic scattering bound- 
ary conditions to the three-component state vector 
gives the S'-matrix elements Sij{e,l, B, A, I). If e > A, 
all three channels are open: i,j = g, p, or e. When 
e < A, as is normally the case in our model, channel e 
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is closed, and Sij is only defined for i,j = g or p. We 
choose the light intensity / low enough that the resuhs 
are in the weak field limit where the Ppg = |S'pgp matrix 
element scales linearly in /. Our results are normalized 
to a standard intensity of / = 1 mW/cm^. 

We find that we can make a change in the asymptotic 
shape of the artificial probe potential to make the model 
much more manageable computationally. The deep po- 
tential well of the excited state and the large kinetic en- 
ergy in the probe channel require a small stepsize in u 
{XAu w 0.005 ao). However, with the large range of u 
(Aumax ~1500-3000 ao) a small Au increases the com- 
putation time and may compromise the numerical sta- 
bility. Therefore, we modify the probe state potential to 
bring Vp{u) to a small negative value at intermediate and 
asymptotic u. This results in a small asymptotic momen- 
tum hk in the probe channel, and allows us to gradually 
increase the stepsize to AAm « 0.5 ao as u increases. The 
coupling between the excited and probe states is turned 
off exponentially before the change in Vp{u). The prob- 
abilities IS'pgp and Ppe are completely independent of 
the asymptotic properties of the probe potential if there 
are no asymptotic barriers. Since we have no centrifu- 
gal potential in the asymptotic probe channel, there are 
no asymptotic centrifugal barriers. The presence of such 
barriers in our previously published model ||20[] resulted 
in some errors at larger I" which we have now corrected 
in the present model. 

The thermally averaged loss rate coefficient via state e 

is: 

K{A,T)^^ r^e-^"^^^ (22) 
hQr Jo ksT 

X + l)\St>g{e,l",B,A,I)\\ 

where Qt = {27r fikBT/h^)^^^ is the translational parti- 
tion function. Identical particle exchange symmetry en- 
sures that only even partial waves exist for the ground 
state. We also define a non-averaged rate coefficient for a 
fixed collision energy £, where we define = s/kB- Only 
the sum over partial waves and branches is performed: 

K{A,e) = ^ J2 {2l" + l)\Sp,{eJ'\B,A,I)\\ 

(23) 

There are two possible cutoffs l^^^y, to the partial wave 
sum provided by the ground and excited-state centrifu- 
gal potentials, respectively. For the ground state we can 
take Z," to be the largest integer for which ?i^^max('max+ 
l)/2/ii?^ < e at the Condon point Re- Thus, the Con- 
don point is classically accessible for /" < I'J^-^^^ and clas- 
sically forbidden for > ^J^ax- For the excited state, 
the centrifugal potential may create a barrier inside the 
Condon point for the g e excitation. The position and 
the height of the barrier depend on J'. For collision en- 
ergies around e = ksTo this barrier may prevent allowed 



ground-state partial waves from contributing to the loss, 
because the excited-state population never reaches the in- 
ner zone where RE decay and SC take place. In this case, 
l'^^^ may decrease from the value defined by the above 
inequality. We find that the ground state cutoff applies 
except for the case of high energy and small detuning. 
In either case \Spg\'^ decreases many orders of magnitude 
as /" varies from l'^-^^^ over the next few /"-values. The 
upper limit for the sum in Eqs. ( [2^ ) and ( p3| ) is set to the 
/"-value where (2/"-h l)|S'pg(e, /", A, /)pis lO'^ of the 
maximum previous {21" -\- l)\Spg{e, I", B, A, /)p value. 

B. Factorization of trap loss probability 

The factorization in Eq. (|^) allows us to separate the 
physics of the long-range excitation and the short-range 
decay to the trap loss channel. Reference shows how 
to determine the short-range probability Ppe from a dif- 
ferent coupled channels calculation where the complex 
decay term —ir/2 in Eq. ( pl| ) is omitted. Since Ppe is 
determined in a region near Rp where the local kinetic 
energy is very high in relation to e, this probability is 
nearly independent of e over a wide range. Therefore, we 
calculate |2^ 

PpeiJ') = \Speis > A, /", B = 0, A, / 0)|2. (24) 

Here e is taken above the threshold energy A where the 
e channel becomes open and Spe is defined. 

Our numerical calculations show, as expected, that Ppe 
is independent of e over a wide range, typically of e/ks 
from 0.3 mK to 300 mK at low J' and 3 mK to 300 
mK at high J', and also independent of A in our small 
range of detuning. The Landau-Zener interpretation of 
Ppe in Eq. (H) leads us to expect that Ppe will vary with 
J'. This variation should be stronger for the outer ^H^- 
■^Sg crossing than for the inner ^'Eu-^Hu crossing. For 
the latter crossing, our calculations do give Ppe values 
which vary slowly with J'. We calculate Ppe{J' = 1) 
to be 0.024, 0.44, 0.31, 0.34 and 0.64 for Mg, Ca, Sr, 
Ba, and Yb respectively. These probabilities are all large 
(order unity) except for the case of Mg. This qualitative 
conclusion is likely to be robust, even though our model 
calculations are only quite approximate. 

In contrast to our results for the ^Etj-'^n„ cross- 
ing. Fig. ^ shows that the calculated Ppe{J') values for 
the outer ^Hg-^T,g crossing indeed depend much more 
strongly on J'. A test of the Landau-Zener formula for 
Mg shows that the result of Eq. (^ is indistinguishable 
from the calculated line on the figure. The qualitative 
feature of a dip in Ppe as it goes near zero for some J' 
is associated with the phase factor in the LZ formula, 
Eq. (H). Since the specific J'-range where this dip occurs 
is sensitive to the potentials used [^, our model calcu- 
lations can only be a qualitative guide even for Mg. The 
relative values for the other species are also only qualita- 
tive guides, since other curve crossings are also involved. 
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In any case, Sr is likely to have a large (order unity), 
perhaps the largest, Ppe for the ^Ilg SC process. 




100 120 



FIG. 4. Calculated probabilities Ppe{J') versus J' for the 
crossine. 

We can use Ppg and Ppe from the close coupling calcu- 
lations to divide out the inner zone probability so as to 
define a numerical excitation-transfer function from the 
long-range region |2|] 



Jeg(£,r,B,A,/) 



Ppg{eJ",B,A,I) 

PpeiJ') 



(25) 



where J' — I" + B. J^g may be interpreted as the prob- 
ability of reaching the short-range region due to opti- 
cal excitation at long-range and propagation to short- 
range, including return after multiple vibrations across 
the short-range well. This interpretation follows from 
the fact that one gets the total trap loss probability 
Ppg{e, I", B, A, /) by muhiplying J^gie, I", B, A, /) by the 
probability Ppe{J') in Eq. ( p4|) of a trap loss event in a 
single complete cycle across the well. 

Using Eq. (^5|), we can define an excitation-transfer 
rate coefficient K{A,e) 



l)Jge{e,l",B,A,I). (26) 



This rate coefficient k(A, e ) is related to the ordinary rate 
coefficient K{A, e) in Eq. (23) through a mean inner zone 
probability {Ppe{e)), which we can define by the relation 



KiA,e) 



{Ppe{s))K{A,e) 



(27) 



Clearly, we can also define a thermal average k{A,T) 
analogous to that in Eq. (^2|), and define a thermal aver- 
age (Ppe(T)) = K{A, T)/k{A, T). 

The usefulness of the factorization in Eq. (|^) is that it 
allows us to define an excitation-transfer rate coefficient 
K from which the inner zone SC probability has been re- 
moved (however, see the discussion in Section III C 3 be- 
low about how a large Ppe{J') may affect the width of res- 
onance features) . We can predict much more confidently 
the properties of the long-range excitation and vibration 



than we can the short-range SC probabilities. Thus, once 
we have a better knowledge of these short-range proba- 
bilities, either through measurements or through better 
theoretical knowledge of potential curves and couplings, 
we can multiply our k coefficients by {Ppe{T)) to get the 
SC rate coefficients. 



C. Limiting cases of the excitation-transfer 
probability 

The attractive molecular potentials support molecu- 
lar vibrational levels with vibrational quantum numbers 
V, as described in Section 



II C . We can find simple 



analytic expressions for the excitation-transfer function 
Jeg = JeePeg factored according to Eq. (||) for two limit- 
ing cases: (1) strongly overlapping resonances where the 
probability is large for spontaneous decay during a sin- 
gle vibrational cycle, i.e., the level width is larger than 
the level spacing, and (2) non-overlapping, or isolated, 
resonances, where many vibrations occur during a vibra- 
tional decay lifetime, i.e., the level width is much smaller 
than the level spacing. For Group II species, tran- 
sitions at small detuning tend to be of the former type, 
but never become fully isolated in the detuning range we 
consider. On the other hand, -'^Ilg transitions tend to be 
of the latter type unless the detuning is very small or the 
SC probability is very large. 



1. Small detuning and fast spontaneous decay 

The quantum mechanical theory of the first limiting 
case for trap loss for small detuning and fast radiative 
decay has been worked out in detail in Refs. |21 4l] , |4^ , p^ , 
where: 



Jeg(e,r',B, A,/) = Jee{e,l'\B,A)Peg{eJ",B,A,T). 



The factor 



Peg{eJ",B,A,I) = 1 



-27rA 



A 



\Veg{Rc)\VinvcDc), 



(28) 



(29) 



where vc and Dc are the speed and slope difference at 
the Condon point, represents the Landau-Zener proba- 
bility of excitation from the ground state to the excited 
state in a one-way passage through the Condon point 
at Re- In this limit, radiative decay is faster than the 
vibrational period {Ty 3> hv^), there are no multiple vi- 
brations, and the excited state transfer factor 



Jeeie,l",B,A) « 1, 



aout = ^ du- 



(30) 
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represents the probability of survival along the clas- 
sical trajectory from the Condon point of excitation 
to the point Rp of inner zone curve crossing; v{u) = 

1 /2 

{2[e — Ve{u) — A]//x} ' is the local classical speed. 

Note that the Landau-Zener expression for P^g in 
Eq. ( p9| ) does not have the proper Wigner law threshold 
behavior, since Jeg should be proportional to k at low col- 
lision energy. However, our numerical Jeg will have the 
proper Wigner law form. Note also that the quantum 
mechanical calculations in Ref. pl| , p4| , ^l|j4^ ] support this 
semiclassical picture of localized excitation at the Con- 
don point, not the delocalized excitation picture of the 
Gallagher-Pritchard (GP) model which for small de- 
tuning predicts a dominant contribution to trap loss from 
off- resonant excitation at distances much less than Re- 
The GP model also does not satisfy the Wigner law at 
low T. We defer detailed comparisons with semiclassical 
theories to a future publication. 



8. Non-overlapping resonances 

The second limiting case is that of non-overlapping vi- 
brational resonances, that is, the spacing hi/^ [see Eq. (j^)] 
between vibrational levels v is much larger than their to- 
tal width Fy. This is typical of large detuning. Then 
jiS'pgp is given by an isolated Breit- Wigner resonance 
scattering formula for photoassociation lines |p§;H : 



r r 

^ vp^ vg 



pa 



[£-(s„ + s,)]' + (r„/2) 



2 ■ 



(31) 



Here Ey — A — Sy is the detuning-dependent position of 
the vibrational level in the molecule-field picture relative 
to the ground state separated atom energy (when A = 
Ey, then Ey — and the vibrational level is in exact 
resonance with colliding atoms with zero kinetic energy) , 
Sy is a level shift due to the laser-induced coupling, and 
the total width = F. 



vp 



■ vg 



v,rad 



is the sum of 

the decay widths into the probe (F„p) and ground state 
(Tyg) channels and the radiative decay rate (Ty^rad)- In 
the weak decay limit (F^, ^ hi/y ), w e can write the Fermi 
golden rule decay widths as |45|]4^ 



F„ = 2TT\{v\Vy^\e,l)\^ = hl^yPy., 



(32) 



where i = g or p, I is the partial wave for channel i, and 
Pyi represents the probability of decay during a single 
cycle of vibration from level v to channel i. For the SC 
process, Pyp is a very weak function of energy as long 
as the detuning is not too large, and we can take Pyp — 
Ppe, where Ppe i s the energy-independent SC probability 
discussed in Sec. HI B . In the weak- field limit, Tyg is very 
small in relation to Ty^rad, and we can ignore it (that is, 
there is no power broadening). 

Using Eq. ( ^2| ) in Eq. (pT|), we get the factorization in 
Eq. (Eq) with the resonant-enhanced transfer function 



J ei 



(hVy) 



[e-iEy+Sy)f + iry/2f' 



(33) 



If we use the reflection approximation for the Franck- 
Condon factor in F„g (34|j3|Q, then 



Peg = 47r^\Veg{Rc)\^ — \M£,l",Rc)\ 



(34) 



Equation ( |34| ) may be used throughout the whole cold 
collision domain (mK to nK). It satisfies the the correct 
Wigner threshold law behavior at low energies because of 
the factor. In the s-wave limit for low temperature, 
we may take the asymptotic form of the ground state 
wavefunction and obtain: 



P.. 



{Rc)? 



hv^Dc 



sin^ k{Rc - Aq) 



(35) 



This looks just like the Landau-Zener result in Eq. (||), 
except that the asymptotic speed Voo appears in the de- 
nominator instead of vc [ p4[ , and the correct quantum 
phase appears in the sine factor instead of a semiclassical 
phase. 



3. Contributions to the linewidths 



The expression, Eq. (|28|), for J^g in the limit of small 
detuning and fast decay does not depend in any way on 
the short-range SC probability. However, in the expres- 
sion for Jeg in the isolated resonance limit, the width F„ 
in the Jee factor, Eq. (|3^), does depend on Ppe through 
the contribution of F^p — TiVyPpe. As long as F^p is small 
compared to Ty^rad, the total width F^, is determined pri- 
marily by Ty^rad^ and the shape of trap loss spectral lines 
will still be nearly independent of Ppe- However, if F„p 
makes a significant contribution to the total width, the 
long-range excitation-transfer function Jeg will show ad- 
ditional broadening dependent on the magnitude of Ppe- 

The total radiative decay width Ty^yad can be calcu- 
lated from the long-range form of the decay rates in 
Eqs. ffl), using the excellent semiclassical approxima- 
tion |4^, (u|F(u)|u) = Uy §^{T{u) /v{u))Xdu, where the 
semiclassical integral is over a complete vibrational cy- 
cle. When Rc < A, we can use the lead term in the 
expansion of T{u) in u in Eqs. (^, so that 



'^v,radC^u) = 2Fat = constaut. 



(36) 



^^*20^''^ 



0.701Fc(^H3), (37) 



where as is defined after Eq. (H) and Fc(^Hg) is evalu- 
ated at the outer turning point of the vibration, which is 
almost the same as the Condon point. 

For the detuning range we consider, the radiative width 
of ^S]„ levels, 2Fat, is much larger than Typ, which can 
be calc ulated from Eq. ( |3^ ) using the probabilities listed 
in Sec. IIIB. Thus, F^ « F„_rad so that the shape of 
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features (that is, their spacings and widths) should be 
well-determined in our calculations. 

Figure H shows Ty^rad and Typ for ^Ilg features for Mg, 
Ca, and Sr. In our detuning range, Ty^rad ^ ^vp for Mg. 
Thus, the shape of Mg ^Hg features should also be well- 
determined in our calculations. On the other hand, for 
Ca and Sr, the Typ is larger due to the larger Ppe. Typ 
increases as A^^^ due to the Vy factor in Eq. (^2|), and be- 
comes larger than Ty^rad near A — SFat in our model for 
Sr and near 20rat for Ca. Thus, we can expect predisso- 
ciation broadening of ^Ilg features to become observable 
for Ca or Sr at relatively small detunings. Measurements 
of such widths could lead to experimental information 
about Ppe for the ^Ilg state. On the other hand, our 
calculated model line shapes should only be viewed as a 
qualitative guide in a region where Typ ^ Ty^rad- 
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FIG. 5. Radiative width Ty^rad for Mg and widths r„p 
Mg, Ca and Sr versus A for the ^Ilg-^Eg SC crossing. 



for 



culating the thermally averaged rate. The main contribu- 
tion to Pdocay comes from the long-range region where the 
potential is determined by its analytic long-range form. 

The Pdocay probability is insensitive to collision energy 
and detuning. In the detuning range A/Fat from 1 to 50 
and for a collision energy of ksTo, we find that Pdocay 
ranges from 0.157 to 0.144 for Mg, 0.103 to 0.100 for Ca, 
0.147 to 0.142 for Sr, 0.113 to 0.110 for Ba, and 0.149 to 
0.145 for Yb. These hardly change at all at a collision 
energy of fcsTc/lOOO, for example, changing to 0.105 to 
0.100 for Ca. 

We have also used a calculation with an artificial probe 
state crossing the excited state potential at Rp {Rp « 
150 an fo r our Mg model for J' = 1), as described in 



Sec. [I G , to calculate the excitation-transfer function J, 



eg 

appropriate to the RE process. We find, as expected, 
that the numerical Jeg function calculated this way is 
very nearly the same as the one calculated using the SC 
Rp at much shorter range. For our detuning range the 
radiative contribution to the total width F^, of ^ E„ levels 
is much larger than contribution due to predissociation 
to the SC channel. 

We expect that our radiative escape trap loss calcu- 
lations are reliable in magnitude, since only long-range 
properties are relevant in determining both J^g and 
^'docay Therefore, in the next Section we can confidently 
give absolute magnitudes for the RE contribution to the 
total trap loss rate coefficient K{A, T) for all alkaline 
earth atoms we study here. 



IV. RESULTS 



A. Trap loss for Mg at T = Td 



D. Radiative escape calculations 

The calculation of the RE trap loss rate coefficient fol- 
lows the factorization procedure in Eq. (Q) as for the 
SC process. The RE loss is not due to a single curve 
crossing, but rather to excited state emission from the 
distance range R < Rp = Xup, where Up is the point for 
which a kinetic-energy increase of 1 K for the atom pair 
has been gained after excitation (the 1 K is arbitrary- 
we only choose it to represent a "standard" loss energy) . 
Clearly, RE can only be significant for the ^T^y state, be- 
cause of the negligible short-range emission from the ^Ilg 
state. We calculate the total probability of radiative es- 
cape, Ppe — ^dccay during a complete cycle of vibration 
across the region u < Uphy integrating along the classical 
trajectory: 

^'dccay(e, J', A) = 1 - exp(-a), 

a = 2A r" du^^. (38) 

-Pdocay depends only weakly on A, J' . Variations with e 
at the highest collision energies also play a role when cal- 



Our calculated results for To = 1.9 mK for Mg colli- 
sions are shown in Figs. |6|(a), 0(a), and ^(a). These re- 
sults are different from the results presented in Ref. , 
since we have corrected some errors we made in the sum 
over partial waves in that reference [Q. Figure |^(a) 
shows on a logarithmic scale the separate contributions 
of each SC or RE process to the thermally averaged rate 
constant K{A, T) from Eq. (|2^), whereas Fig. |7|(a) shows 
the corresponding results for A'(A, e) at a single collision 
energy s. Figure ^(a) shows on a linear scale the sum 
of contributions from all loss processes, and shows what 
one might expect to see in a laboratory spectrum. 

The dominant loss process for Mg at 1.9 mK is due to 
RE from the -^E„ state. The spectra for RE and SC 
processes have the same shape, since they have the same 
excitation-transfer function n. The RE and SC processes 
differ only by a multiplicative factor that is nearly in- 
dependent of A, due to the different Ppe factors for RE 
and SC. The RE probability only varies by 0.157 to 
0.144 from detunings of 1 to 50 Fat, whereas the SC 
probability is constant over this range. The spec- 
tra in Figs. |6|(a) and M(a) are nearly the same, since the 
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FIG. 6. Contributions from the ^Eu RE and ^Eu and ^Tlg 
SC processes to the thermally averaged loss rate coefficient 
K(A,T) at (a) 1.9 mK and (b) 190 pK as a function of 
laser detuning A for Mg at a standard laser intensity 1 = 1 
mW/cm^. The scales for the excitation-transfer coefficients, 
K(A,r), for the SC processes are indicated by the vertical 
axes to the right. The vibrational quantum numbers from 
the top of the potential, vt, are indicated for ^Eu and ^Ilg 
features. 




FIG. 7. Contributions from the ^E^ RE and ^E„ and ^Ug 
SC processes to the loss rate coefficients K{A, e) at a fixed 
collision energy (a) e = fcs(1.9 mK) (b) e = ^^(190 /iK) as a 
function of laser detuning A for Mg at laser intensity 1 = 1 
mW/cm^ K{A,e) is a sum over partial waves and branches 
for e = fcflT. The corresponding excitation- transfer coeffi- 
cients, k(A, e), are indicated by the vertical axis to the right. 
Excited state rotational quantum numbers are indicated for 
the Vt = 24 ^Ilg feature in (a). 



broad features do not change much upon thermal averag- 
ing. The rate coefficient becomes very small as detuning 
decreases below 2 or 3 Fat- This is because Jee «C 1 for 
very small detuning due to spontaneous emission during 
the long-range approach of the two atoms. Spontaneous 
emission losses become small for detunings larger than 
around 10 Fat, and vibrationally resolved, but rotation- 
ally unresolved, photoassociation structure begins to de- 
velop as detuning increases. This occurs as the spacing 
between adjacent vibrational levels from Eq. be- 
comes larger than the radiative decay width. Several ro- 
tational features with different J' may contribute to each 
of the broad photoassociation resonances, with the range 
of J' depending on detuning. Each individual rota- 
tional line has a width on the order of (2Fat + fcsT)//i w 
200 MHz. There is negligible predissociation broadening 
due to SC processes in thi s reg ion of the spectrum. 

Using Eq. (|^) in Section [I C , the vibrational quantum 
number vt , as counted down from the top of the potential 
at the dissociation limit, can be given for the resolved, 
or partially resolved, features in the trap loss spectrum. 
We define vt to he vjj — v rounded up to the next inte- 



ger. Each integer vt value defines an energy range which 
contains only one vibrational level for a given J. The 
Vt quantum numbers for ^S„ and ^Hg features are indi- 
cated on Fig. ^. Note that there are many levels (not 
calculated) within the range A/Fat < 1, a range where 
Eq. is not meaningful due to retardation effects on 
the potential. The broad ^E„ features provide an exam- 
ple of overlapping resonances, analogous to those treated 
by Bell and Seaton |Q for the case of dielectronic recom- 
bination where the spacing between collisional resonance 
levels becomes less than their radiative decay width. 

The contribution to K{A,T) from the ^Hg SC pro- 
cess shows much sharper vibrational structure than the 
corresponding ^S„ spectrum. This is because of the 
small radiative widths of the ^Ilg levels, which become 
even smaller as A increases. The individual contribution 
from a number of narrow rotational levels is evident in 
Fig. 0(a). Figure ||(a) shows that this ^Ug structure even 
survives thermal averaging. Figure ||(a) shows that sharp 
^Ilg features can even survive thermal averaging at 1.9 
mK, although such features are quite weak for Mg and 
would be hard to see (However, see below for Ca and 
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FIG. 8. Total thermally averaged Mg spectrum, K{A, T) 
summed over all RE and SC contributions, on a linear scale, 
(a) At Td =1.9 mK, (b) In the s-wave limit at T = 1.9 /xK. 
The vibrational quantum numbers Vt are indicated for the 
and ^Ilg features. Only excited J' — 1 levels contribute 
R-branch transition from s- waves in panel (b). 



Ba, where such features might be observable). We find 
that there are sharp J' = 1 ^Ilg features due to s-wave 
collisions that can be much narrower than ksT (which is 
about 40 MHz at 1.9 mK), whereas features due to /" > 
collisions have widths on the order of fc^T. This s-wave 
behavior is evident in our numerical calculations, but can 
be easily explained in terms of the analytic behavior of 
the isolated fine shapes using Eqs. (^, (|3^), and (34). 
We will discuss this s-wave resonance narrowing feature 
in more detail elsewhere. 

Figures || and ^ both show that at very small detuning, 
on the order of 1 or 2 Fat, the trap loss is dominated 



by SC due to the ^Hg state. The increasing radiative 
transition probability as detuning decreases, and the near 
absence of spontaneous emission losses for the weakly 
emitting state, ensures that the ^Ilg contribution to trap 
loss must be dominant at very small A. We will show 
in the next section that this is even more important for 
the heavier species. Our conclusion concerning the role 
of the ^EL state at small A agrees with the findings of 
Refs. [0j2gL 

Figures 3|(b) and 0(b) show the contributions to SC 
and RE processes for Mg at 190 /iK. The broad fea- 
tures are not very sensitive to changing the temperature. 
They narrow slightly at the lower temperature. However, 
the ^Hg features simplify and clearly have contributions 
from fewer partial waves. The effect of thermal averag- 
ing on ^Hg features is to cause some broadening, with 
consequent decrease in peak height. 



<i 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 

: (a) 




1 , 


1 1 1 1 1 1 1 1 

Mg 'n 

T= 1.9 ^IK 






ir , 

r v.rad 











10.35 10.4 10.45 10.5 10.55 10.6 
A/r , 



H 
<i 




10.35 



10.6 



10.4 10.45 10.5 10.55 

A/r , 
at 

FIG. 9. Single ^Ilg vibrational feature in the vicinity of 
A lOFat for (a) Mg, (b) Ca, and (c) Sr. The figure shows 
the quality of the isolated resonance approximation for the 
excitation-transfer line shape K(A,r). The solid line is the 
complex potential numerical calculation, and the dashed line 
is the analytic line shape based on Eqs. (0), (0), and (^sj). 



B. Trap loss for Mg near 1 /iK 

Figure ^(b) shows if(A, T) summed over all compo- 
nents at the extremely cold temperature of 1.9 /iK. This 
is deeply in the Wigner law domain, where only s-wave 
collisions contribute to the spectrum, and the rate con- 
stant -fC(A, T) becomes independent of T [|l[. The broad 
features are similar to the ones at higher tempera- 
ture, but are due only to absorption by a single R branch 
line from Z" = to a J' = 1 level. The only signifi- 
cant broadening is due to radiative decay. On the other 
hand, the ^ 

line from /" = to a J' 
nent sharp features in the spectrum, having widths on 
the order of a few MHz due to radiative decay. Even the 



Hg features, also due to a single R branch 
1 ^Hg level, become promi- 
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level near 1 Fat detuning is quite sharp and isolated. Sec- 



tion HE discusses why we expect to get the vibrational 
spacings right, although we do not expect to predict cor- 
rectly the actual position of levels, which depend on an 
unknown phase due to the short-range ^Ilg potential. 

Figure ^a) shows the excellent quality of the isolated 
resonance approximation for a Mg ^FEg s-wave absorp- 
tion feature due to a single vibrational level. This ap- 
proximation should be good in this case, since the mean 
vibrational spacing near this level is 280 MHz, which is 
much larger than the width. The Figure compares the 
numerical line shape with that calculated using the iso- 
lated resonance formulas discussed in Section |[II C| . The 
analytic formula calculates the factors in Eqs. (p8D, which 
are used in Eq. ( Pw, b y making the isolated resonance ap- 
proximation, Eq. (|33|), and the reflection approximation, 
Eq. (|35|). The linewidth in the denominator of Eq. (33), 
calculated to be 1.6 MHz from Eq. (p7|), is almost entirely 
due to weak spontaneous decay of this ^Hg level, as dis- 
cussed Section [H C in relation to Fig. ^. Any broaden- 
ing due to thermal averaging is negligible, since ksT/h ~ 
0.04 MHz. 

We have verified that our thermal spectrum at rela- 
tively high temperature, 1.9 mK, is to a good approxima- 
tion independent of the choice of ground state potential, 
as discussed in Section II D. This is because of the need 



to sum over several partial waves, for which the phase 
of the ground state wavefunction varies by more than tt. 
In addition, the need to average over a range of collision 
energies also contributes a range of phase variation to the 
ground state wavefunction. 

The spectrum at very low temperature, on the other 
hand, is sensitive to the phase of the ground state wave- 
function, which is generally unknown for Group II species 
and strongly dependent of the details of the ground state 
potential. This sensitivity is explained by the reflection 
approximation in Eq. (|35[), which shows Peg is propor- 
tional to siia^ k{Rc — Aq). We have just seen that the 
reflection approximation is excellent for isolated reso- 
nance line shapes. Therefore, if we know K{A,T) in 
the s-wave domain for one scattering length Aq, and 
if we have a different potential with a different scat- 
tering length the K{A,T) for the new case can be 
scaled from the original one by multiplying by the ratio 
sin^ k{Rc — Aq)/ sin^ k{Rc — Aq). Figure |o| compares 
this scaling (dashed lines) to numerical calculations (solid 
lines) for several different model ground state potentials 
with different Aq. The former are scaled from our original 
calculation, for which Aq = —95 ap. Figure |l^ demon- 
strates that this scaling is a good approximation, even 
when the scattering length is unusually large and even for 
overlapping features. The scaling relation is excel- 
lent at small A for scattering lengths havin g magnitudes 
up to a few times xq (defined in Sec. II D and having a 
value of 36 ao for Mg) . The scaling is even a reasonable 
approximation for the case where Aq = 400 ao and the 
ground state wavefunction has a node at Rc = Aq near 



A/Fat = 15. The node for the Aq = 930 ao case occurs 
for A/Fat < 1 and is off scale in Fig. |l^ for the A'q = 99 
ao case. 
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FIG. 10. Scaling with different scattering lengths of 
K{A,e) for the transition in Mg. The bold solid line 
shows the numerically calculated K{A,e) at e = fcfl(1.9 
fiK) for the "standard" ground state model potential with 
Aq = —95 ao. The other solid lines show the calculated 
K{A, e) for three other model potentials with different scat- 
tering lengths of 99 ao, 400 ao, and 930 ao. The dashed line 
shows the scaled K(A, e) calculated from the "standard" one 
using the scaling relation discussed in the text. The detuning 
for which the Condon point is 400 ao is indicated by the ar- 
row. The effect of the node in the ground state wavefunction 
is evident for the Aq — 400 ao case. 



V. OTHER ALKALINE EARTH ATOMS 

Our calculations for the other Group II species and Yb 
are shown in Figs. |(b), |(c), |ll|, |l|, and ]l3[ We trust 
that these model calculations, which can only provide 
order of magnitude estimates for SC probabilities and 
predissociation contributions to linewidths, will provide 
a useful qualitative guide to differences and similarities 
among the various cases to guide future experiments on 
these systems. Our calculations should be fairly robust 
with respect to qualitative expectations as to the different 
kinds of features to expect in trap loss spectra. 

Figure |^(b) shows that a very low temperature Ca 
^Hg feature is very similar to the Mg one previously 
discussed. The total width is slightly larger than the 
radiative width due to weak predissociation of this fea- 
ture (see Fig. ||) . The effect of the large predissociation 
width, where F„„ > Frad, is evident for the Sr feature in 



- vp 

Fig. |9|(c). The isolated resonance approximation is also 
beginning to fail for Sr lines because of the strong predis- 
sociation broadening in our model with Ppe = 0.22 (see 
Fig. |4|) . Although our model calculation for Sr predissoci- 
ation widths should not be considered to be reliable, the 
model does show that if the widths of ^Hg features like 
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FIG. 11. Excitation transfer coefficients K(A,e) on a loga- 
rithmic scale for the (a) ^'Eu and (b) ^Ilg states as a function 
of laser detuning A for Mg, Ca, Sr, Ba, and Yb at a laser 
intensity of 7 = 1 mW/cm'^. 



the one in Fig. ^(c) could be measured, the data should 
allow a value to be determined for Ppg. Since temper- 
atures in the nK regime have already been reported for 
intercombination line cooling of Sr, it may be quite fea- 
sible to measure such widths. 



Figure 11 shows the thermally averaged excitation- 
transfer coefficients k(A,T) (see Eq. ( p6| ) and following) 
for the and ^Ilg states in these systems at To for the 
^S-^ cooling transition. In spite of the fact that the 
inner zone SC probability Ppe is divided out of the expres- 
sion for k(A, r), there are still a number of differences 
among the different species. The differences in spacing 
and contrast of the individual vibrational features that 
appear at larger detuning is clearly related to the vibra- 
tional spacing, Eq. (^), which decreases with increasing 
mass. The differences in magnitude can be qualitatively 
related to the scaling of the different factors that make up 
k(A, e) in Eq. (Eq). There are four factors that contribute 
to the scaling: (1) 1/Qt fJ.~^^^, (2) the sum over 
I" - ^Lx - Mdo^VA^/^ (3) Jee(peak) (z.jr.)2 ^ 
X'A'/y{^^dl'/^), and (4) P^g ^ \V,giRc)\V (vDc) ^ 
/i^/^dp^'^/A^^'^- The net scaling of the peak magnitude 
of k{A,T) thus scales approximately as X^/{pdQ^^A^/'^). 
This gives scaling factors at the same A of 1, 5.0, 3.4, 
6.4, 1.1 for Mg, Ca, Sr, Ba, and Yb respectively (these 



factors should be scaled by an additional factor of X/d^^^ 
if evaluated at the same scaled detuning, A/Tat)- These 
scaling factors account for the relative magnitudes of the 
peak k{A,T) for the ^S]„ state in Fig. |l^(a) in the rel- 
atively flat region from 20 to 50 A/Fat. The scaling for 
the ^Ilg spectra in Fig. |ll|(b) also needs to take into ac- 
count the predissociation contribution to the width Ty, 
which was taken to be purely radiative for the scaling of 
the spectrum in Fig. |l^(a). For example, this extra 
predissociation broadening lowers the peak of Sr features 
below those for Mg in Fig. |ll|(b). 
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-Ca 830|iK[10''2cm^/sl 
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FIG. 12. Spectrum K{A, T) summed over all RE and SC 
contributions at Td for Mg, Ca, Sr, and Ba at a laser intensity 
of / = 1 mW/cm^ 

Figure ^ shows our model thermally averaged 
K{A, T) summed over all contributions. With the caveat 
that the relative contributions of SC processes are not 
likely to be reliable in our model calculations, these model 
spectra show qualitative features that one might observe 
in laboratory spectra. In particular, ^Ilg SC processes 
make a dominant contribution to the small detuning trap 
loss for A < a few Fat- This has already been discussed 
in Refs. We can compare our results to the mea- 

sured 2K{A,T) = 4.5(0.3)(1.1) x 10'^° cm^/s |§ for 
Sr at A/Fat = 1.75, / = 60 mW/cm^, and T « 4Td 
[Q. Although the effect of a strong laser field needs 
to be investigated for this case, Ref. suggests that 
near-linear scaling may apply to small detuning trap loss 
even in the strong field domain (see Fig. 6 of that refer- 
ence). If we assume linear scaling with /, our calculated 
value for T — Tjj at / = 1 mW/cm^ scales to a value of 
2K = 6x 10"^° cm^/s at / = 60 mW/cm^. The agree- 
ment of our very approximate model to within a factor 
of two with the measured result for Sr is gratifying and 
fends confidence to the usefulness of our estimates. 

At the present, there are no other data on Sr or other 
Group II species to which we can compare our calcula- 
tions directly. The Ca2 photoassociation spectra in a 3 
mK Ca MOT reported by Zinner et al. in Ref. ex- 
tend over a detuning range from about 50 to 2700 Fat, 
which is larger than we calculate. They observed ^E^ 
features and gave a detailed analysis of partially resolved 
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rotational substructure for a feature near A = 27 GHz 
= 780 Fat- The fact that the width of this feature could 
be explained by a combination of natural and thermal 
broadening of several rotational lines implies that pre- 
dissociation broadening makes a small contribution to 
the lincwidth of this feature. If we assume that 20 MHz 
or less of the observed 150 MHz feature width is due 
to predissociation, we would then estimate the SC 
Ppe < 0.05, which is much less than the value 0.4 esti- 
mated by our model for Ca. Although our model should 
not be extrapolated to such large detuning without care- 
ful testing, this apparent inconsistency points out that 
much more detailed knowledge of potentials and matrix 
elements is needed for accurate calculations. It is an in- 
teresting fact to be explained why the apparent predisso- 
ciation rate of levels in Ca2 is relatively small, given 
the likelihood of several curve crossings with moderately 
large matrix elements (see Fig. ||). 

Our calculations in Fig. |lj suggest that resolved struc- 
ture due to ^Hg features may be seen at small detunings 
below around 25Fat for Ca and Ba. Structure for Sr 
is predicted to be suppressed by strong predissociation 
broadening. No ^H^ structure was reported for detun- 
ings larger than around 50Fat in Ref. Such ^Hg 
structure in Ca.2 at these larger detunings may be hard 
to see due to masking by the strong features. 

Figure |l^ shows our predictions for Ca and Sr features 
at extremely low T = T/j/lOOO. This is in the s-wave 
limit where the ^Hg structure becomes quite sharp, as 
discussed in relation to Fig. ^ above. In this domain 
sharp ^Hg features should be the dominant features in the 
trap loss spectrum. It is noteworthy that this structure 
is predicted to persist even to very small detunings on 
the order of Fat. Thus, if the Sr trap loss experiments 
of Ref. [0 could be repeated at these low temperatures, 
such features might be measurable. Figure ^ predicts 
that predissociation widths may be large enough for Sr2 
"'^Hg features at even a few Fat detuning that observed 
broadening in the spectra might be able to determine 
Ppe for the Sr ^Hg SC process. Thus, low temperature 
measurements provide for tests of consistency with high 
temperature measurements. 

VI. CONCLUSIONS 

We have carried out model calculations of the small- 
detuning coUisional trap loss spectrum of laser cooled 
Group II species Mg, Ca, Sr, and Ba and also Yb. We 
consider detunings A up to 50 atomic linewidths Fat to 
the red of the ^Sq ^Pi laser cooling transition for 
these species and treat both inelastic state-changing col- 
lisions and radiative loss. Although our calculations are 
only model calculations because the short-range molec- 
ular potentials are not known to sufhcient accuracy, we 
do incorporate the correct long-range aspects of the po- 
tentials and spectra. These calculations are intended as 
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FIG. 13. Spectrum K{A, T) summed over all RE and SC 
contributions at s-wave domain at Td/IOOO for (a) Ca and 
(b) Sr at a laser intensity of 7 = 1 mW/cm'^. 

a guide for developing experimental studies on these sys- 
tems, which have the advantage that the collisions are 
not complicated by molecular hypcrfine structure. 

We consider both the mK range for Doppler cooling 
on the allowed ^So — ^ ^Pi transition, and /iK range for 
Doppler cooling on the ^Sq — > '^Pi intercombination tran- 
sition. Collisions in the mK range involve many partial 
waves, whereas /iK collisions only involve s-wave colli- 
sions. Our quantum mechanical calculations avoid semi- 
classical approximations and properly account for the 
threshold properties of the collisions. Our interpretation 
of trap loss collision dynamics is based on a factorization 
of the overall probability into parts that represent long- 
range excitation, propagation to the short-range region, 
and short-range radiative or curve crossing processes that 
lead to loss. Thus, we can define an excitation-transfer 
coefficient k(A,T), which, unlike the conventional rate 
coefficient K{A, T), offers a significant degree of indepen- 
dence from the details of unknown short-range processes. 
Our analysis shows how analytic formulas in the limits of 
small or large detuning can be used to interpret the trap 
loss spectrum. 

The trap loss spectra in all the Group II systems 
are influenced by two molecular transitions, the dipole- 
allowed ^Sg — > ^'Eu transtion and the dipole- forbidden 
^Sg ^Hg transition. The latter becomes allowed at 
long-range because of retardation corrections to the tran- 
sition matrix element. The features are structure- 
less at small detuning and reduced in magnitude due to 
spontaneous decay of the excited state as the atoms ap- 
proach one another on the excited state molecular poten- 
tial. They show broad vibrationally resolved but rota- 
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tionally unresolved photoassociation structure as detun- 
ing increases away from atomic resonance. On the other 
hand, the ^Ilg absorption always dominates at small de- 
tuning. Resolved ^Ilg vibrational and rotational pho- 
toassociation structure can persist even to small detun- 
ing, and should be especially prominent at very low tem- 
perature. Measurement of the widths of such features 
could lead to information about the short-range prob- 
ability of the state-changing collisions, at least for the 
heavier Group II elements. 

There are only very limited data with which we can 
compare our calculations. Our model calculations agree 
within a factor of two with the measured Sr trap loss 
rate coefficient at a single detuning. Photoassociation 
spectra for Ca only exist for much larger detuning than 
we consider here, but suggest that the probability of 
the state changing process may be much smaller 
than our model calculations indicate. The time is right 
for more detailed and complete experimental studies on 
these Group II systems. Recent experimental advances 
in Group II cooling and trapping suggest that such stud- 
ies will be forthcoming. A number of other directions 
are also open for continuing experimental and theoret- 
ical studies, for example, trap loss collisions near the 
^So ^Pi intercombination line, or colhsions associated 
with two ^Pi atoms or two ■^P atoms. 
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